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ABSTRACT 



In this paper, we present a multi-color photometric study of the nearby spiral galaxy M81, 
using images obtained with the Beijing Astronomical Observatory 60/90 cm Schmidt Telescope 
in 13 intermediate-band filters from 3800 to 10000A. The observations cover the whole area of 
M81 with a total integration of 51 hours from February 1995 to February 1997. This provides 
a multi-color map of M81 in pixels of 1" .7 x 1".7. Using theoretical stellar population synthesis 
models, we demonstrate that some BATC colors and color indices can be used to disentangle 
the age and metallicity effect. We compare in detail the observed properties of M81 with the 
predictions from population synthesis models and quantify the relative chemical abundance, age 
and reddening distributions for different components of M81. We find that the metallicity of 
M81 is about Z = 0.03 with no significant difference over the whole galaxy. In contrast, an age 
gradient is found between stellar populations of the central regions and of the bulge and disk 
regions of M81: the stellar population in its central regions is older than 8 Gyr while the disk 
stars are considerably younger, ~ 2 Gyr. We also give the reddening distribution in M81. Some 
dust lanes are found in the galaxy bulge region and the reddening in the outer disk is higher than 
that in the central regions. 



Subject headings: galaxies: abundances - galaxies: evolution - galaxies: individual (M81) - 
galaxies: stellar content - (ISM:) dust, reddening 
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INTRODUCTION 



Spatially resolved information about the age, metallicity and interstellar medium reddening of galaxies 
is a powerful tool to study galaxy evolution since it provides essential clues in star formation history chemical 



composition and enrichment history and environment of galaxies (Buzzoni 1989). To obtain such information 



we need to know the stellar component and the overall properties of stellar populations ( Lcithcrcr & Hcckman 
19951 ) . Ideally one would like to study resolved individual stars in galaxies. However, given the limited spatial 



resolution of current telescopes, this is only possible for a few very nearby galaxies. As a result, the stellar 



content of even some relatively simple galaxies remain to be unraveled (Thuan 1991). Information of stellar 



population and star formation history in these galaxies, however, can still be obtained from studying the 



integrated properties of the stars (e.g. Bchmitt et al. 1996, Goerdt & Kollatschny 1998) 



Since the pioneering work of Tinsley (1972) and Searle et al. (1973), evolutionary population synthesis 
has become a standard technique to study the stellar populations of galaxies. This is a result of improvements 
in the theory of the chemical evolution of galaxies, star formation, stellar evolution and atmospheres, and 
the development of synthesis algorithms and the availability of various evolutionary synthesis models. A 
comprehensive compilation of such models was published by Leitherer et al. (1996) and Kennicutt (1998). 
Widely used models include those from the Padova and Geneva group (e.g. Schaerer fe dc Kotcr 199*^ , 



Schaerer fe Vacca 1998, Bressan et al. 1996, and Chiosi et al. 199Sj), GISSEL96 (Chariot fe Bruzual 



1991| , Pruzual fe Chariot 1993| , |Bruzual fc Chariot 1996; ), PEGASE ( |Fioc fc Rocca-Volmerange 1997p and 



STARBURST99 ( [Leitherer et al. 1999|) 



Many previous studies of integrated stellar populations use spectroscopic data, usually for limited regions 
in galaxies. In this paper, we will, instead, use multi-color photometry to probe the stellar populations. The 
multi-color photometry provides accurate spectral energy distributions (SEDs) for the whole galaxy, although 
at low spectral resolution. We shall demonstrate that it is a powerful tool to study the structure and evolution 
of the galaxy together with the theoretical evolutionary population synthesis methods (for an application of 
a similar technique with, but with fewer colors, to moderate redshifts, see Abraham et al. 1999). For this 
purpose, we pick M81 as the first application of this multi-color approach. 

M81 is an excellent candidate because it is a nearby early-type Sab spiral galaxy at a distance of 3.6 
Mpc and with an angular size of ~ 26'. The angular extent is large enough such that the disk and bulge 
regions are well separated from the ground. It has been the subject of numerous previous studies providing 
a wealth of information with which to compare the new metallicity and internal reddening distribution. The 
internal reddening has been studied by Kaufman et al. (1987, 1989), Devereux et al. (1995), Ho et al. (1996) 
and Allen et al. (1997). The metallicity has been studied by Stauffer & Bothun (1984), Garnett & Shields 
(1987), Brodie & Huchra (1991) and Perelmuter et al. (1995). In this paper, we present a further detailed 
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study of M81 using the unique dataset obtained from the BATC [] multi-color sky survey. 

The outline of the paper is as follows. Details of observations and data reduction are given in section 
2. In section 3, we provide a brief description of the model, and analyze the evolution of the integrated 
colors, color indices with age and metallicity. The observed two-dimensional spectral energy distributions 
(SEDs) of M81 were analyzed using stellar population synthesis models of Bruzual & Chariot (1996). The 
distributions of metallicity, age and interstellar reddening are given in section 4. In section 5, we discuss how 
different star formation histories and stellar population synthesis models change our results, and compare 
our results with previous studies. . The conclusions are summarized in section 6. 



2. OBSERVATIONS AND DATA REDUCTION 

2.1. CCD Image Observation 

The large field multi-color observations of the spiral galaxy M81 were obtained in the BATC photometric 
system. The telescope used is the 60/90 cm f/3 Schmidt Telescope of Beijing Astronomical Observatory 
(BAO), located at the Xinglong station. A Ford Aerospace 2048x2048 CCD camera with 15/zm pixel size is 
mounted at the Schmidt focus of the telescope. The field of view of the CCD is 58' x 58' with a pixel scale 
of 1".7. 

The multi-color BATC filter system includes 15 intermediate-band filters, covering the total optical 
wavelength range from 3000 to 10000A (see Fan et al, 1996). The filters were specifically designed to avoid 
contamination from the brightest and most variable night sky emission lines. A full description of the 
BAO Schmidt telescope, CCD, data-taking system, and definition of the BATC filter systems are detailed 
elsewhere ( Fan et al. 1996] , Chen ct al. 200C). To study the age, metallicity and interstellar reddening of 



M81, the images of M81 covering most part of the optical body of M81 were accumulated in 13 intermediate 
band filters with a total exposure time of about 51 hours from February 5, 1995 to February 19, 1997. The 
CCD images are centered at RA = 09 h 55 m 35 s .25 and DEC=69°21'50".9 (J2000). The dome flat-field images 
were taken by using a diffuse plate in front of correcting plate of the Schmidt telescope. For flux calibration, 
the Oke-Gunn primary flux standard stars HD19445, HD84937, BD+262606 and BD+174708 were observed 
during photometric nights. The parameters of the filters and the statistics of the observations are given in 
Table 1. 



^The Beijing- Arizona- Taiwan-Connecticut Multicolor Sky Survey 
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Table 1: Parameters of the BATC filters and statistics of observations 
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"Central wavelength for each BATC filter 
b Image numbers for each BATC filter 

c Zero point error, in magnitude, for each filter as obtained from the standard stars 
^Background errors before image smoothing 
e Background errors after image smoothing 

2.2. Image data reduction 

The data were reduced with standard procedures, including bias subtraction and flat-fielding of the CCD 
images, with an automatic data reduction software named PIPELINE 1 developed for the BATC multi-color 
sky survey (Fan et al. 1996). The flat-fielded images of each color were combined by integer pixel shifting. 
The cosmic rays and bad pixels were corrected by comparison of multiple images during combination. The 
images were re-centered and position calibrated using the HST Guide Star Catalogue. The sky background 
of the images was obtained by fitting image areas free of stars and galaxies using the method described in 
Zheng et al. (1999). The absolute flux of intermediate-band filter images was calibrated using observations 
of standard stars. Fluxes as observed through the BATC filters for the Oke-Gunn stars were derived by 



convolving the SEDs of these stars with the measured BATC filter transmission functions (Fan et al. 1996). 
Column 6 in Table 1 gives the zero point error, in magnitude, for the standard stars in each filter. The 
formal errors we obtain for these stars in the 13 BATC filters is <J 0.02 mag. This indicates that we can 
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define the standard BATC system to an accuracy of <J 0.02 mag. 

After background subtraction, the standard deviation of the background for each image is 2.0ADU. 
Because the signal-to-noise ratio decreases from the center to the edge of the galaxy, we smoothed the 
images with a boxcar filter. The window sizes of the boxcar were selected depending on the ADU values 
of the BATC10 band image (7010A). If the ADU value was less than 100, the pixel was set to zero; if the 
value was higher than 100, the pixel was adaptive-smoothed by boxcar filter of N x N (cell size), where 



N=min(151/ VADUbatcio, 15). By this method, the images were smoothed depending on the S/N of each 
cell. In the central area of M81, the original pixels were used, whereas near the edge of M81 the mean value 
of multiple pixels (cells) were used, as a result, the spatial resolution decreased from center to outer edge. 
The background errors before and after smoothing are given in the last two columns in Table 1. 

Finally, the flux derived at each point of M81 is listed in Table 2 Q . The results provide a two-dimensional 
spectral energy distributions (SED) for M81. The ADU number of each image were converted into units 
of 10 _30 ergs _1 cm -2 Hz -1 . As an example, we present some SEDs for the different areas of M81 in this 
paper. The table gives the following information: Column 1 and Column 2 give the (X, Y) positions of the 
photometric center of the regions, in units of arcseconds. The coordinate system is centered on the nucleus of 
the galaxy (RA = 09 h 55 m 35 s 25; DEC = 69°21'50".9, in J2000). The A-axis is along the E-W direction with 
positive values towards the east, and the Y~-axis is along the N-S direction with positive values towards the 
north. Column 3 to Column 14 give the fluxes relative to the BATC08 filter (6075A). Column 15 gives the 
flux in the BATC08 filter in units of 10 _30 ergs _1 cm~ 2 Hz~ 1 . For convenience in later discussions, we define 
the 'central' areas as regions with |ARA| < l'.O, |ADEC| < 1'.9, the 'bulge' region with 1.0 < |ARA| < 
3'.4, 1.9 < |ADEC| < 4'.0 and the 'disk' region with 3'.4 < |ARA| < 6'.2,4'.0 < |ADEC| < 8'.0. 

We show a black and white image of M81 in Figure 1. To display the features clearer, a "true-color" image 
of M81 is available in the electronic version, which is combined with the "blue" (3894A), "green" (5785A) and 
"red"(7010A) filters. The filters selected here are free from any strong emission lines. From this image, we 
can see directly the stellar population difference in different areas of the spiral galaxy M81. In the following 
sections, we will analyze quantitatively the stellar populations in M81 with our 13 color data. 

The bright HII regions consist of young clusters, and the evolutionary population synthesis methods 
used in this paper do not represent young clusters well. In addition, the central nucleus of M81 exhibits some 
of the same characteristics as classical Seyfert galaxies, has no evidence of stellar clusters or a population 
of hot young stars ( Kaufman ct al. 1996 , Devereux et al. 1997 , Davidge fc Courteau 1999| ). So we mask 
some bright HII regions and the central nucleus of M81 (shown as white spots in Fig. 4a, 5, 6, including the 
foreground stars) in this paper. 



2 The full Table 2 and color versions of Figs. 1, 4a, 5 and 6 are available electronically. 
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Tablc 2. Two dimensional spectral energy distributions for the central, bulge and disk areas of M81 
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Note. — Rows 1-5, 6-10, 11-15 are for the central, bulge and disk areas respectively. The absolute flux (last column) 
is in units of 10 _30 ergs _1 cm _2 Hz _1 . The full table for M81 is available electronically. The pixel values for the center of 
M81 is (300, 375). 




Fig. 1— "True-color" estimate of the M81 generated by using the BATC02 (3894A) filter image for blue, 
BATC07 (5785A) for green, and BATC10 (7010A) for red; the image is balanced by making the background 
old population orange and hot stars blue. The center (origin) of the image is located at a = 09 h 55 m 35 s .25, 5 = 
69°21'50".9(J2000.0). The image size is about 17'.0 by 21'.0. North is up and east is to the left. We refer to 
the region with |ARA| < l'.O, |ADEC| < 1'.9 as the central region, the region with 1.0 < |ARA| < 3'.4, 1.9 < 
|ADEC| < 4'.0 as the bulge, and the region with 3'.4 < |ARA| < 6'.2,4'.0 < |ADEC| < 8'.0 as the disk. A 
black and white image of Figure 1 be displayed in paper, the "true-color" image of Figure 1 is available in 
the electronic version. 
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3. DATABASES OF SIMPLE STELLAR POPULATIONS 

A simple stellar populations (SSP) is denned as a single generation of coeval stars with fixed parameters 



such as metallicity, initial mass function, etc (Buzzoni 1997). In evolution synthesis models, they are mod- 



eled by a collection of stellar evolutionary tracks with different masses and initial chemical compositions, 
supplemented with a library of stellar spectra for stars at different evolutionary stages. Because SSPs are the 
basic building blocks of synthetic spectra of galaxies that can be used to infer the formation and subsequent 



evolution of the parent galaxies ( Jablonka et al. 1996 ). In order to study the integrated properties of stel 



lar population in M81, as the first step, we use the SSPs of Galaxy Isochrone Synthesis Spectra Evolution 
Library (Bruzual & Chariot 1996 hereafter GSSP). We study the SSPs as the first step for two reasons. 
First, they are simple and reasonably well understood, so it is important to see what one can learn using 
this simplest assumption, and then check whether more complex star formation history give qualitatively 
similar conclusions. This is a common approach often taken in the evolutionary population synthesis models 



for galaxies ( Vazdekis et al. 1997 , Mayya 1995 ). Second, although we assume each pixel is described by an 
SSP, we emphasize that the whole galaxy is not SSP; so our assumption is not as strong as it may seem. 
Nevertheless, this is a significant assumption. Fortunately, it appears the adoption of more complex star 
formation history does not change the results qualititatively; we return to this issue in §5.1. 



3.1. Spectral Energy Distribution of GSSPs 

The Bruzual & Chariot (1996) study has extended the Bruzual & Chariot (1993) evolutionary population 
synthesis models. The updated version provides the evolution of the spectrophotometric properties for a 
wide range of stellar metallicity. They are based on the stellar evolution tracks computed by Bressan et al. 
(1993), Fagotto et al. (1994), and by Girardi et al. (1996), who use the radiative opacities of Iglcsias et 
al. (1992). This library includes tracks for stars with metallicities Z = 0.0004,0.004,0.008,0.02,0.05, and 
0.1, with the helium abundance given by Y = 2.5Z + 0.23 (The solar metallicity is Zq — 0.02). The stellar 
spectra library are from Lejeune et al. (1997,1998) for all the metallicities listed above, which in turn consist 
of Kurucz (1995) spectra for the hotter stars (O-K), Bessell et al. (1991) and Fluks et al. (1994) spectra 
for M giants, and Allard & Hauschildt (1995) spectra for M dwarfs. The initial mass function is assumed to 
follow the Salpeter's (1955) form, dN/dM oc M -2 - 35 , with a lower cutoff M\ = 0.1M Q and an upper cutoff 



Af u = 125M Q (Sawicki fc Yee 199S 



3.2. Integrated Colors of GSSPs 



To determine the age, metallicity and interstellar medium reddening distribution for M81, we find 
the best match between the observed colors and the predictions of GSSP for each cell of M81. Since the 
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observational data are integrated luminosity, to make comparisons, we first convolve the SED of GSSP with 
BATC filter profiles to obtain the optical and near-infrared integrated luminosity. The integrated luminosity 
L\ { (t, Z) of the zth BATC filter can be calculated with 

t (tz) jwgj^M^ m 
A mln (i) ^( x > dX 

where the F\(t,Z) is the spectral energy distribution of the GSSP of metallicity Z at age t, <Pi(X) is the 
response functions of the BATC filter system, and A m i n (i) and A max (i) are respectively the maximum and 
the minimum effective wavelength of the zth filter (i = 1,2, • • •, 13). 

The absolute luminosity can be obtained if we know the distance to a galaxy and the extinction along 
the line of sight. Since we do not know the exact distance to M81, in this paper, we shall work with the 
colors that are independent of the distance. We calculate the integrated colors of a GSSP relative to the 
BATC filter BATC08 (A = 6075A): 

C Xi {t,Z) =L Xi (t,Z)/L 6075 (t,Z). (2) 

As a result, we obtain intermediate-band colors for 6 metallicities from Z = 0.0004 to Z = 0.1. In the 
panels of Fig. 2, we plot the colors as a function of age for GSSP with different metallicities. The following 
remarks can be made, (a) It is apparent that there is a uniform tendency for SSPs to become redder for all 
colors as the metallicity increases from Z = 0.0004 to Z = 0.05. The near-UV and optical colors show the 
same qualitative behavior as those at longer wavelengths, (b) There is a wide range in age (from 1 to 20 Gyr) 
in which the colors vary monotonically with time except for the highest metallicity Z = 0.1. Therefore, once 
we know the metallicity and interstellar reddening, we can use these colors to determine the age distribution 
of M81, provided that the stellar population is well modeled by SSPs. (c) For Z = 0.1, there is only a 
limited age range for the monotonic behavior in colors. One reason for this behavior is the appearance of 
AGB-manque stars at Z = 0.1. These stars skip the AGB phase and directly go through a long-lived hot 



HB phase (Bruzual fc Chariot 1996). There are very few, if any, examples of galactic stars with such a high 



metallicity. So our results are not affected by this peculiar high-metallicity case. 

3.3. Color Indices of GSSPs 

The observed colors are affected by interstellar reddening, which will of course complicate our inter- 



pretations ( pstlin ct al. 1998| ). The interstellar reddening in the center region of M81 can be measured 



by its emission lines, but for the outer regions, the problem becomes very complex. If we suppose that 
the extinction law from 3800A to 10000A has no high frequency features, the spectral indices will not be 
affected much by the uncertainties in the extinction, so we can use the spectral indices to reduce the effect 
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Age (Gyr) Age (Gyr) Age (Gyr) Age (Gyr) 



Fig. 2. — Evolution of eight selected colors C\ t (cf. eq. for simple stellar population (SSP) models 
for different metallicities as predicted by the GISSEL96 library. Each line pattern represents a different 
metallicity, as indicated inside the frame. All the models shown are computed with a the Salpeter(1955) 
initial mass function with a lower cutoff of 0.1M Q and an upper cutoff of 125M©. 
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of interstellar extinction. 

Spectral indices are (by definition) constructed by means of a central band pass and two pseudo- 
continuum band-pass on either side of the central band. The continuum flux is interpolated between the 



middle-points of the pseudo-continuum band passes (Bressan et al. f996, Worthcy et al. 1997). Since we 



only observed M81 with intermediate-band filters, not a genuine one-dimensional spectrum, so we must use 
some pseudo-color indices to replace the conventional definitions. We define a color index I\ j (t, Z) of a SSP 

by 

h.(t,Z)=L Xj (t,Z)/L x . +1 (t,Z), (3) 

where L\ j (t, Z) is the luminosity of a SSP with the metallicity Z at age t and wavelength Aj, L\ j+1 (t, Z) is 
the luminosity in the (j + 1 )th filter for the same SSP. The color indices can reduce the effect of reddening, 
especially in the wavelength region longer than 5000A. 

Among all the BATC filter bands, we find that the color index centered at 85fOA (Js5io) is much 
more sensitive to the metallicity than to the age; the center of this filter band is near the Call triplets 
(A = 8498, 8542, 8662A). The strength of Call triplet depends on the effective temperature, surface gravity 
and the metallicity for late type stars ( ^hou 1991 ) . 



In an old stellar system, the effect of metallicity on Call triplet becomes prominent. In fact, we find 
that there is a very good relation between the flux ratio of issio = ^8510/^9170 and the metallicity for stellar 
populations older than 1 Gyr. We plot this relation in Figure 3. Similar relations are also found in many 
other observation and stellar population synthesis models. 

The relation shown in Figure 3 is crucial for our metallicity determination and later studies, so it is 
important to check whether this is indeed a reliable method. An early study by Alloin & Bica (1989), based 
on the analysis of stars, star clusters and galaxy nuclei indicated a strong correlation of Call triplet with the 
surface gravity, logg. However, further studies suggested that the Call triplet strength depends not only 
on the surface gravity but also on the metallicity. A detailed analysis of the behavior of the Call triplet 
feature as a function of stellar parameters was performed by Erdelyi-Mendes & Barbuy (1991), making use 
of a large grid of synthetic spectra. They concluded that Call triplet has a weak dependence on the effective 
temperature, a modest dependence on surface gravity, but a quite important dependence on metallicity. 
They even suggested that the Call triplet strength may vary exponentially with the metallicity. Moreover, 
these lines have been studied by Diaz et al. (1989) , Mallik (1994) and Idiart et al. (1997). They have also 
suggested the Call triplet strengths depend on the metallicity (Mayya 1997). So the relation between the 
785io and metallicity seems to be reliable and can be used to determine the metallicities; our own investigation 
of GSSPs seems to be consistent with these recent studies. 
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Fig. 3. — a) Dependence of the pseudo color-indice Tssio with metallicity for the GISSEL96 SSP models. 
The line symbols are the same as in Fig. 2. The mctallicities are labeled above each line, b) The correlation 
between the pseudo color-index Issio an d metallicity, for five different ages (T = 15,10,5,3,1 Gyr). The 
solid line shows the best fit through the curves. The lower panel for the PSSPs, the upper panel for the 
GSSPs. 
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4. DISTRIBUTION OF METALLICITY, AGE AND REDDENING 



In general, the SED of a stellar system depends on age, metallicity and reddening along the line of 



sight. The effects of age, metallicity and reddening are difficult to separate (e.g., Calzetti 1997, Origlia et 



al. 1999, Vazdekis et al. 1997). Older age, higher metallicity or larger reddening all lead to a redder SEDs of 
stellar systems in the optical ( Molla et al. 1997 , Brcssan et al. 1996| ). In order to separate the effects of age, 
metallicity and interstellar reddening of M81, we first determine the metallicity by the color index IssiO) as 
discussed above, and then obtain the age and reddening by using GSSP model of known metallicity and a 
extinction law (see §4.2). 



4.1. Metallicity distribution 

As discussed in §3.3, there is a good correlation between the color index Issio and the metallicity. We 
will use the relation obtained from GSSP; this relation is similar for other stellar population synthesis models. 
We find that the correlation can be fit with a simple formula: 

Z = (0.83 -0.84 x 7 85 io) 2 . (4) 

This curve is shown in the upper panel of Figure 3. The scatter in Figure 3 is due to the difference in age. 
If the ages are younger than 1 Gyr, the scatter becomes larger. Figure 3 shows that for a stellar system of 
ages older than 0.5 Gyr, we can estimate the metallicity with an error less than the interval of metallicity 
given by GSSPs. 

Using this method we obtained the metallicities for each part of M81 except for the nucleus and the Ha 
line emission region (we have masked them out, see §2.2). Figure 4a shows the resulting metallicity map of 
M81. Figure 4b shows the radial distribution of the metallicity, the curve is derived from the Fig. 4a by 
averaging over ellipses of widths 17" along the major axis. We used an inclination angle of i = 59° and a 
position angle of PA = 157° for the major axis of the galaxy. 

To our surprise, we do not find, within our errors, any obvious metallicity gradient from the central 
region to the bulge and disk of M81. In most parts of M81, the mean metallicity is about 0.03 with variation 
<J 0.005. These results are identical to past suggestions that early-type spirals may have relative high 
abundances and weak gradients. Taking into account of the age scatter, the true value of the metallicity is 
likely within a range between Z = 0.02 and Z = 0.05. From the metallicity map of M81, we can also clearly 
see that in some outer regions the metallicities are higher; most of these regions are located in spiral arms 
and around HII regions, where a younger stellar population is present. 
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Fig. 4a. — A map of the population metallicity in M81. Light gray corresponds to the low and dark gray to 
the high metallicity. Some bright HII regions, the central nucleus of M81 and most of foreground stars are 
masked, shown as white spots. 
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Fig. 4b. — The radial distribution of metallicity for M81, averaging over ellipses along the major axis (its 
position angle is taken to be 157°). The line is the regression line obtained by a polynomial fit. No obvious 
metallicity gradient exists from the central region to the bulge and disk of M81. 
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4.2. Age and reddening distribution 

Since we model the stellar populations by SSPs, the observed colors for each cell are determined by two 
parameters: age, t, and dust reddening, E(B — V). In this section, we will determine these parameters for 
M81 simultaneously by a least square method. The procedure is as follows. For given reddening and age 
(recall that the metallicity is known, see the previous subsection), we can obtain the predicted integrated 
colors by convolving the dust-free predictions from GSSP with the extinction curve given by Zombeck (1990). 
The best fit age and reddening values are found by minimizing the difference between the observed colors 
and the predicted values: 

12 

i? 2 (x, y, t, Z, E) = YF&^i V) - Z ' ( 5 ) 

i=l 

where C^ p (t, Z, E) represents integrated color in the ith filter of a SSP with age t, metallicity Z and reddening 
correction E, and Cf* s (x,y) is the observed integrated color at position (x,y). 

Figure 5 shows the age map for M81. It clearly indicates that the stellar population in the central 
regions is much older than that in the outer regions and the youngest components reside in the spiral arms 
of M81. There is a smooth age gradient from the center of the galaxy to the edge of the bulge. The age 
in the innermost central region (within 17") is older than J> 15 Gyr. The age at the more extended central 
region is about 9 Gyr. In the bulge edge area, the age is about 4 Gyr. In contrast, the stellar component in 
the disk area is much younger than that in the bulge region. The mean age in disk area is about 2.0 Gyr. 
We can see that the age in spiral arms is even younger than the inter-arm areas, about 1 Gyr. 

Because the age obtained in the outer disk region is around 1 Gyr, the metallicity, which is determined 
by the color indices, might have large errors (see §4.1). The errors in turn will make the age determination 
uncertain. Therefore the age for disk can only be regarded as a rough estimation. However, the general 
trend in the age distribution should be reliable. 

Figure 6 shows the reddening map for M81. From this figure, we find a large difference in reddening 
between the bulge and the disk. In the bulge, the reddening, E(B-V), is in the range of 0.08 to 0.15. For 
regions where E(B-V) is larger than 0.1, we found some spiral like cirrus. There is a very obvious high 
reddening lane around the nucleus with reddening equal to or higher than the disk area. The maximum 
of E(B-V) in the lane reaches 0.25. This half-loop lane can be verified in the future with IR or CO line 
observations. In the disk area, the mean reddening of E(B-V) is about 0.20. In the central regions, the mean 
reddening of E(B-V) is very small. These results suggest that dust is largely absent in the central regions of 
M81. The dust seems to be distributed mainly along the inner part of spiral arms and around the nucleus. 
Color versions of Figs. 4-6 are available electronically. 
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Fig. 5. — A map of the population ages in M81. Light gray corresponds to the young and dark gray to the 
oldest zones. Same as in Fig. 4a, the white spots represent the masked regions, such as bright HII regions, 
the central nucleus of M81 and most of foreground stars. 
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5. DISCUSSION 

The results we presented so far are based on the (strong) assumption that all stars in a small region form 
in an instantaneous burst and hence the stellar population of each cell can be modeled as SSPs. Unfortunately, 
the star formation rate history is an essential but very uncertain ingredient in the evolutionary population 
synthesis method since it can vary from galaxy to galaxy and from region to region inside a single galaxy 
as well. It is only for simplicity that we have adopted an instantaneous star formation history, clearly it is 
important to check whether the results are significantly changed if one varies the star formation history; we 
address this issue in §5.1. While there seems to be general agreement between the GISSEL96 SSPs models 



(which we used) with other similar ones (Chariot et al. 1996), there are some fine differences. In §5.2, we 



study how the results are changed if we adopt a different population synthesis model from the Padova group. 
In §5.3, we compare our results with earlier works. 

5.1. Continuous Star Formation 

We assume that stars are formed from the interstellar gas exponentially with a characteristic time-scale 



r, i.e., = ^oexp t ' T . This model is often used to calculate the integrated colors of galaxies (Kcnnicutt 



199q) and allows a more diverse star formation history. If r — > oo the model approximates constant star 



formation rate, while for r — > it approximates an instantaneous burst (Abraham et al. 1999). It seems 



that spiral galaxies could be well fitted with r of the order of several Gyr (Fioc et al. 1997). 

Since we do not know the appropriate r value for M81, we have explored the values r = 0.1, 1, 3 Gyr. 
We have calculated the age, metallicity and interstellar reddening distributions for each value of r. As the 
value of r increases, the age (t) increases throughout M81, while the interstellar reddening decreases and 
the metallicity is little changed. Although the numerical values of these quantities do change in each point 
of M81, the two dimensional distributions of age, metallicity and interstellar reddening of M81 for different 
values of r are similar to the ones shown in Figs. 4-6. 



5.2. Comparison With Other SSP Models 

There exist a number of SSP models that are synthesized with different approaches. It is important 
to check the sensitivity of our results to the SSP models adopted. The SSP models from the Padova 
group (hereafter, PSSPs) is suitable for the comparison, because they use a similar technique of "isochrone 
synthesis" to predict the spectral evolution of stellar populations. The PSSPs provide the basis for the 



population synthesis models (Bressan et al. 1994, see Bressan et al. 1996 for revisions and extensions) 



The PSSPs use a comprehensive set of stellar evolutionary tracks of the Padova group for a wide range of 
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initial chemical compositions from Z = 0.0004 to Z = 0.1 with AY/AZ = 2.5. The initial masses of the 
evolutionary tracks cover the range of 0.6— 120M Q , except for the set of metallicity Z=0.1, where the masses 
are from 0.6 — 9M Q . The initial mass function is the Salpeter (1955) law. More details can be found in 
Bressan et al. (1994), Silva (1995), and Tantalo et al. (1996). The main difference between GSSPs and 
PSSPs is the library of stellar spectra: the GSSPs use theoretical stellar spectra from Lejeune et al. (1997) 
while the PSSPs use theoretical stellar spectra from Kurucz (1992). 

Using our method, we calculate the colors and color indices for each PSSP in the BATC filter system. 
Using similar procedures as for the GSSP (see §4), we obtained the metallicity, age and reddening distri- 
butions of M81. For PSSP, the metallicity map of M81 again has no obvious metallicity gradient, but the 
mean metallicity is somewhat higher, about 0.035. There is a smooth age gradient from the center of M81 
to the edge of the bulge, except that the mean age in the disk area is lower than 1 Gyr. The interstellar 
reddening value from PSSP is obvious bigger than that from GSSP. In the bulge, the reddening value is in 
the range of 0.18 to 0.35. In the disk area, the mean reddening of E(B-V) is about 0.40. The reddening 
value in the central region amount to 0.15. The distributions of metallicity, age and interstellar reddening 
are very similar to those found using GSSP. 



5.3. Comparison with previous work 

Numerous determinations of the amount of extinction in M81 have recently been obtained. Allen et 
al. (1997) compared the detailed distribution of HI, Ha, 150 nm Far-UV continuum emission in the spiral 
arms of M81. They found every reliable bright peak in the Ha has a peak in the Far-UV, and concluded 
that the effects of extinction on the morphology are small on the spiral arms. Filippenko & Sargent (1988), 
based on the ratio of the narrow components of Ha and H(3, concluded that the central regions of M81 is 
reddened by E(B-V)=0.094 mag. These results are very similar to our results for internal reddening, that 
the mean reddening in the spiral arms and in the central regions of M81 are small. In addition, Kaufman et 



al. (1987, 1989), using the Ha and radio continuum (Bash & Kaufman 1986) observations, have studied the 



distribution of extinction along the spiral arms in M81, and obtained a mean A v = 1.1±0.4 mag for 42 giant 
HII regions with high surface brightness. Hill et al. (1992) used Near-UV, Far-UV and V-band images of 
M81 and cluster models to derive an A v — 1.5 mag for the HII regions on the arms. These internal reddening 
values are larger than ours for spiral arms. It can be explained for two reasons. First, the internal reddening 
value of Kaufman et al. are derived from the giant HII regions of M81 with high surface brightness, but 
these bright HII regions are excluded in our study. So we have have used different regions for the internal 
reddening study, it is therefore not clear that they must agree. Second, the continuum regions are affected 
less by reddening than the emission line regions in a galaxy. The internal reddening value of Kaufman et 
al. is the reddening for emission line regions, but our results are from the continuum regions. The former 
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is significantly larger than the latter; this systematic difference has been seen in many other emission line 



galaxies ( {Kong fc Cheng 1999| , |Calzetti 1997[ ). 

Perelmuter et al. (1995) have obtained spectra for 25 globular cluster in M81. Following the method of 
Brodic & Huchra (1990), based on the weighted mean of six indices, they measured these clusters' metallicity. 
The mean metallicity was calculated both from the weighted mean of the individual metallicities, and directly 
from the cumulative spectrum of the 25 globulars. Both results yielded the same value, (Fe/H) = — 1.48±0.19 
(Z = 0.033), which is identical to that derived by Brodie & Huchra (1991) using 8 clusters ((Fe/H) = 
— 1.46 ± 0.31). No correlation has been observed between magnitude and metallicity of the globulars in the 
Milky Way and in M31. Thus the mean metallicity of the 25 globulars should be representative of the M81 
system as a whole. These results agree with our results for the metallicity of M81 (Z « 0.03) very well. 
On the other hand, using the low-dispersion spectra of 10 HII regions in M81, Stauffer & Bothun (1984) 
estimated the oxygen abundances for those HII regions from the observed emission lines. They derived a 
mean abundance near solar and a weak abundance gradient in M81. Using the empirical calibration method 
and the photoionization models, Garnett & Shields (1987) have analyzed the metallicity abundance and 
abundance gradients for 18 HII regions in the galaxy M81. The major result of this study is the presence 
of order-of-magnitude gradient in the oxygen abundance across the disk of M81. These results differ from 
ours. However, the differences can be readily explained. First, as for the internal reddening, the previous 
results are derived from the bright HII regions of M81, the abundance gradient is therefore for these HII 
regions. But our result come from the whole galaxy except these bright HII regions. So it is not clear that 
they should have the same behaviors in metallicity. Another caveat to the previous results is that the Pager 
et al. (1980) abundance calibration (which was used in Stauffer & Bothun 1984, Garnett & Shields 1987) is 
not expected to be exact. 

At last, we must emphasize that although the method we used in this paper can be used to constrain the 
variation of metallicity, population age, and reddening across M81 are for the central region, the bulge, and 
the disk minus the spiral arms, it may not be suited to study the property of the spiral arms. There are two 
main reasons. First, there are hundreds of HII regions on the spiral arms, and, the evolutionary population 
synthesis methods that we used in this paper are not well represented very young clusters. Second, the 
signal-to-noise ratio decreases from the center to the edge of the galaxy, we have smoothed the images with a 
boxcar filter (See §2.2). For the outer disk regions, such as the spiral arms, the smoothing tends to blend the 
hundreds of HII regions with their surroundings. We will study these bright HII regions at the spiral arms 
of M81 in more detail in a subsequent paper when we have better data, using the evolutionary population 
synthesis method, and show how well it can work for young clusters. 
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6. Conclusions 

In this paper, we have, for the first time, obtained a two-dimensional SED of M81 in 13 intermediate 
colors with the BAO 60/90 cm Schmidt telescope. Below, we summarize our main conclusions. 

• Using the new extensive grid of GSSPs covering a wide range of metallicity and age, we calculated 
the colors and color indices for 13 colors in BATC intermediate-band filter system. We find that 
some of them can be used to break the age and metallicity degeneracy, which enables us to obtain 
two-dimensional maps of metallicity, interstellar reddening and age of M81. 

• From the two dimensional metallicity distribution of M81, we find no obvious metallicity gradient from 
the central regions to the outer disk. In most part of M81, the mean metallicity is about 0.03 with 
variation <J 0.005. Some regions in M81, however, have higher metallicity; they are mostly located in 
the spiral arms and around HII regions, where the younger component resides. 

• From the two dimensional age distribution in M81, we find that the mean ages of the stellar populations 
in the central regions are older than those in the outer regions, which suggests that star formations in 
the central regions occurred earlier than the outer regions. 

• We find a strong difference in reddening between the bulge region and the disk region. In the bulge 
area, the reddening, E(B-V), is in the range of 0.08 to 0.15. The mean reddening in the disk area is 
higher, about 0.2. There are some high reddening spiral-like cirrus in the bulge. 

• In order to understand how sensitive our method is to different assumptions about star formation 
history and different stellar population synthesis models, we have studied an exponential star formation 
history and compared the results obtained with GSSP and PSSP. We find that although the precise 
values of age, metallicity and interstellar reddening are different, the general trend of in the metallicity, 
age and reddening distributions is similar. 

• Finally, we have compared the internal reddening and metallicity maps of M81 with previous studies. 
We find that the agreements are generally good. In addition, we find that the properties for the bright 
HII regions and other parts may be different. 

The results of M81 presented here illustrate that our method and observational data provide an efficient 
way to study the distribution of metallicity, age and interstellar reddening for nearby face-on galaxies. Similar 
data have already been collected for similar galaxies M13, NGC589 and NGC5055. The analysis results of 
these galaxies will be published in a forthcoming paper. 
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